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ABSTRACT 

We identify two new tidally distorted white dwarfs (WDs), SDSS J174140.49+652638.7 and 
J211921.96-b01825.8 (hereafter J1741 and J2119). Both stars are extremely low-mass (ELM, 
< Q.2Mq) WDs in short-period, detached binary systems. High-speed photometric observations 
obtained at the McDonald Observatory reveal ellipsoidal variations and Doppler beaming in both sys- 
tems; J1741, with a minimum companion mass of 1.1 Mq, has one of the strongest Doppler beaming 
signals ever observed in a binary system (0.59±0.06% amplitude). We use the observed ellipsoidal 
variations to constrain the radius of each WD. For J1741, the star's radius must exceed 0.074 Rq. For 
J2119, the radius exceeds 0.10 Rq. These indirect radius measurements are comparable to the radius 
measurements for the bloated WD companions to A-stars found by the Kepler spacecraft, and they 
constitute some of the largest radii inferred for any WD. Surprisingly, J1741 also appears to show a 
0.23±0.06% reflection effect, and we discuss possible sources for this excess heating. Both J1741 and 
J2119 are strong gravitational wave sources, and the time-of-minimum of the ellipsoidal variations can 
be used to detect the orbital period decay. This may be possible on a timescale of a decade or less. 
Subject headings: binaries: close — Galaxy: stellar content — stars: individual (SDSS 

J174140.49-h652638.7, SDSS J211921.96-001825.8) — stars: neutron — stars: 

variables: general — white dwarfs 



1. INTRODUCTION 

White dwarfs (WDs) are the remnant, degenerate cores 
of stars. About half of the stars in the galaxy go through 
stellar evolution as single stars and end up as typical 
Earth-sized 0.6 Mq C/O-core WDs. The remaining 
half evolve in binary systems, including short-period sys- 
tems where stable or unstable mass transfer can take 
place in the late stages of their stellar evolution. Many 
short-period systems go through one or two common- 
envelope phases and evolve into systems containing low- 
mass WDs. Indeed, radial velocity surveys of low-mass, 
He-core WDs { M < 0.45Mf:i) indicate that most form i n 
binary systems (jMarsh et al.lll99"5t iBrown et al.ll20lTal ). 
The lower mass WDs are predicted to be larger in radius 
than their more massive counterparts. However, obser- 
vational data for radius measurements for low-mass WDs 
are scarce. 

Eclipsing binaries are important laboratories for con- 
straining the mass-radius relation for WDs, but only 
a handful of direct radius measurements for low-mass. 
He-core WDs exist. Observations using the Kepler 
spacecraft have so far found four A-stars with bloated 
~ 0.2 — O.4M0 WD companions. The radius mea- 
surements foi__th£se_Jou^_WDsran£e^ 
0.15 Ra^ (lva,n KerkwiikeFaLl lMot iCarter et al.] [20T1I: 
iBreton et al.l 1201 1[ ). However, due to their extreme 
environments, all four of these WDs are predicted to 
be larger than normal. Hence, they provide only 



an upper limit for the mass-radius relation for low- 
mass WDs. In addition, there have been four eclips- 
ing low-mass WD sy stems detected froi n the ground 
( Stei nfadt et al.ll20ia iParsons et all 1201 It iVeiines et al.l 
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I20rit IBrown et al.ll2011bD . IMaxted et al.l poill ) recentlv 
identified a 0.23Mo stripped core of a red giant star, a 
pre-He WD. They measure a radius of 0.33 Rq from the 
eclipse observations, providing an estimate for the radii 
of the extremely low mass (ELM) WD progenitors. 

Tidally distorted WDs in (non-) eclipsing systems also 
provide reliable constraints on WD radii. The ampli- 
tude of the ellipsoidal variations is roughly S/elv ~ 
{m2/mi){ri/a)^, where a is the orbita l semi-major axis 
and r i is the radius of the primary fe.g-. lMazeh fc Faigleii 
|2010[ ). Hence, this effect is most important for short- 
period systems with extreme mass ratios and more im- 
portantly with relatively large (radius) WDs, i.e., sys- 
tems containing ELM WDs. Not surprisingly, the first 
and second direct detection of tidally distorted WDs oc- 
curred in t he P < 1 h orbit al period E LM WD sys- 
tems J 0106 ()Kilic et al.ll20llbl ) and J0651 ()Brown et al.l 
l2011bf) . J0651, a 12.75-minute orbital period detached 
eclipsing double WD system, illustrates the wealth of 
photometric information that can be present in an ELM 
WD system: there are primary and secondary eclipses, 
ellipsoidal variations, and a strong Doppler beaming sig- 
nal. 

The ELM Surv ev (IBrown et al.l 120101 : iKihc et al.l 
[20TTat IBrown et al.ll2012h continues to discover a large 
number of short-period systems containing low-mass 
WDs. This survey has been quite successful, so far 
uncovering two dozen d ouble-degenerate binaries with 
< 10 Gyr merger times (|Kilic et al.ll2012D and the only 
known detached binary WDs with <1 h orbital peri- 
ods: t wo 39-minute orbit al period binaries, J0106 and 
J1630 (jKilic et al.ll2011bl l3). and the 12.75-minute sys- 
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tern, J0651 (|Brown et al.ll20irbl ). Based on the proba- 
bility of eclipses and other photometric variability at or- 
bital period timescales, we have established a follow-up 
program to observe the shortest-period binaries found in 
the ELM Survey. These time-series photometric obser- 
vations have been carried out at the McDonald Observa- 
tory using the Argos instrument, a frame-transfer CCD 
mount ed at the prime focus of the 2.1m Otto Struve tele- 
scope ()Nather fc Mukadamll200l . 

Observations of tidal distortions in these short-period 
systems can provide important constraints on the phys- 
ical parameters of the binary. Any mis-alignment be- 
tween the WD spin period and the binary orbital pe- 
riod can be measured using the ellipsoidal variations; a 
mis- align ment will result in tidal heat ing of the ELM 
WD (e.g.. |Pirol[20TltlFuller fc LaillMTI ). If the radius of 
the WD is known, ellipsoidal variations can also be used 
to constrain the mass of the unseen companion and the 
inclination angle. This has strong implications for un- 
derstanding the future evolution of these sy stems, which 
depen ds on the mass ratio of the two stars ([Marsh et al.l 

l2onl . 

Here we take advantage of the latest discoveries from 
the ELM Survey. J1741 and J2119 were recently identi- 
fied as 1.5-2 hr orbital period systems containing two of 
the lowes t surface grayity WDs cu rrently known (log g — 
5.2-5.4: [B^own et al.|[20loll2012[ ). We report the detec- 
tion of ellipsoidal variations and Doppler beaming sig- 
nals in both stars. We discuss our use of the Doppler 
beaming signal to confirm the radial velocity variations 
previously observed in J 1741, and use the ellipsoidal vari- 
ations to constrain the radii of the tidally distorted ELM 
WDs. We discuss our observations and analysis of J1741 
and J2119 in Sections 2 and 3, respectively. We summa- 
rize our results and the mass-radius relation for low-mass 
WDs in Section 4. 

2. SDSS J174140.49+652638.7 

J 1741 is one of th e lowest surface gravity WDs known 
()Brown et al.l [20T2| ). It has Teff = 9,790 ± 240 K, 
\ogg = 5.19 ± 0.06, and M - 0.16 M©. J1741 has a 
1.4666 ± 0.0001 h orbital period and a mass function 
of 0.830 ± 0.018 Mq. Given the spectroscopically de- 
termined mass of the primary, the minimum mass of the 
unseen companion is 1.09 Mq, and there is a 57% chance 
that the inclination is such that its companion is more 
massive than 1.4 Mq. 

2.1. The Light Curve 

Photometric observations of J1741 were carried out at 
the McDonald Observatory over two nights in 2011 May 
and one night in 2011 September, for a total of more than 
9.5 hr of coverage. Exposure times for this g — 18.4 mag 
WD ranged from 15 to 20 s, depending on the seeing and 
sky transparency. Observations were obtained through a 
1mm BG40 filter to reduce sky noise. 

We performed weighted aperture photometry on the 
calibrated frames using the external IRAF package 
ccd_hsp written by Antonio Kanaan. We divided the 
sky-subtracted light curves by at least two brighter com- 
parison stars in the field to allow for fluctuations in 
seeing and cloud cover. Usi ng the wqed software suite 
(|ThomDSon fc Mullallvl[2009l) . we fit a second-order poly- 
nomial to the data to remove the long-term trend caused 
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Fig. 1. — High-speed photometry of J1741 over 9.5 hr (top panel), 
folded at the orbital period, 87.998 minutes. The bottom panel 
shows the same light curve binned by 100 phase bins. The solid line 
includes 1.46% amplitude ellipsoidal variations, a 0.59% Doppler 
beaming signal, and a 0.23% reflection effect. 

by atmospheric extinction, and applied a timing correc- 
tion to each observation to account for the motion of the 
Earth around the barycenter of the solar system. 

Figure [1] shows all 2,156 Argos light curve points ob- 
tained for J1741 (top panel) folded over the best-fit or- 
bital period from the radial velocity variations, 87.998 
minutes. The bottom panel shows the same light curve 
binned into 100 orbital phase bins. The red solid line 
corresponds to the best fit orbital solution, as discussed 
in the following analysis. 

2.2. Analysis and Discussion 

There are four major effects that can cause photo- 
metric variability in the primary of a binary system: 
Doppler beaming, reprocessed light from the secondary 
(reflection), ellipsoidal variations, and eclipses. Here, 
and throughout, we will consider the ELM WD to be 
the primary, as it is the only visible component of this 
single-lined spectroscopic binary. 

Given a high enough radial velocity variation, Doppler 
beaming (also referred to as Doppler boosting or rela- 
tivistic beaming) will act to m odulate stellar flux upon 
approach or recession (see, e.g.. lZucker et al.ll2007D . We 
treat this as a sin0 modulation at the orbital period, 
where = is defined by the spectroscopic conjunction 
(the point at which the primary is farthest away from 
us). Given our blue-bandpass filter, a Teff = 9,790 ± 240 
K, and a radial velocity amplitude K = 508 ± 4 km 
s~^, w e expect a 0.58 ±0.0 2% Doppler beaming signal in 
J1741 (IShporer et al.ll201Cil ). When present, this signal is 
usually evident as a strong asymmetry in the maxima of 
ellipsoidal variations. 

Irradiation of the primary by the secondary can cause a 
reflection effect, which will also sample the orbital period, 
with a maximum at = 0. Thus, we treat the reflection 
effect as a cos (j) modulation at the orbital period. 

Finally, tidal distortions of the primary will cause el- 
lipsoidal variations. The dominant modulation occurs 
when the larger face comes into view twice per rota- 
tion, effectively a cos 20 modulation of the spin period, 
which would be the orbital period for a synchronized sys- 
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tern. However, since tidal distortions do not cause a per- 
fectly ellipsoidal shape, we treat the ellipsoidal variations 
as harmonics to the first four cos terms as derived in 
IMorris fc Naftilanl (fl993h . 

To best characterize the variability in the light curve, 
we perform a Monte Carlo analysis. We create 10^ syn- 
thetic light curves by replacing the measured flux / with 
f + g Sf, where Sf is the error in flux and g is a Gaus- 
sian deviate with zero mean and unit variance. We then 
fit each light curve with a five-parameter model that 
includes an offset, and the (co)sine terms for Doppler 
beaming, ellipsoidal variations, refl ection, and the first 
harmonic of the orbital period (e.g.. iShporer et al.|[201Cll : 
ISirko fc Pac^vKskil [2001 IMazeh fc Faiglei< 120101) . Our 
high-speed photometric observations reveal significant 
amplitudes to the sin 0, cos 0, and cos 20 terms (the rest 
are consistent with zero within the errors; our full results 
can be found in Table [T]) . 

Additionally, we have computed a Fourier transform of 
our time-series photometry. With less than 10 hr of data 
spread over four months, there is much alias structure 
around our highest peak, but the third-highest alias is 
identical to half the radial velocity period, within the 
formal error estimates: Pelv = 2,639.920±0.016 s. This 
peak occurs with 1.43 ± 0.06% amplitude. (These are 
formal least-squares errors; if we did not have a guess for 
the period a priori, the spread in alias peaks would make 
a more realistic error, yielding Pelv — 2,667 ± 27 s.) 
After pre-whitening by this first periodicity, the second- 
highest alias of the remaining peak occurs at Pdb = 
5,280.27 ± 0.13 s (0.67 ± 0.06% amplitude). 

However, since we have signals that may occur at cos (j) 
and sin (/), our Monte Carlo analysis yields a more robust 
estimate for the amplitude of the Doppler beaming effect: 
0.59 ± 0.06% amplitude. The Doppler beaming signal 
yields no new physical information about the binary, but 
since it matches the expected value of 0.58±0.02% within 
our errors, it confirms the radial velocity amplitude ob- 
served in this system, as well as the spectroscopically 
determined temperature for the ELM WD. This is also 
by far the strongest measurement of Doppler beaming 
in a binary system from the ground, even greater than 
the 12 .75-minute double WD binary J0651 ([Brown et al.l 
[2011bD , due to its c ooler effective temper ature (see Equa- 
tions (2) and (3) in IShporer et al.|[20100 . 

The strongest modulation measured in our Monte 
Carlo analysis occurs for the cos 2(/) term, consistent with 
ellipsoidal variations, with an amplitude of 1.46 ±0.06%. 
These variations are extremely useful for constraining the 
radius of the tidally distorted ELM WD. To do so, we use 
the discrete Fourier series for ellipsoidal variations from 
[Morris & Naftilan (1993). Equation (1) in that work 
yields a theoretical prediction for the ellipsoidal varia- 
tion amplitude dominated by 



L(0)/Lo 



~3(15 + iii)(l -|-Ti)(ri/a)3gsin^i 
20(3 -Ml) 



cos(20) 



based on the mass function q ~ m2/mi, the limb- 
darkening (ui) and gravity-darkening (ri) coefficients for 
the primary, the semimajor axis of the system a, the or- 
bital inclination i, and the radius of the primary ri. Us- 
ing reasonable values for the limb-darkening and gravity- 
darkening coefficients, we constrain the radius of the pri- 
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Fig. 2. — Allowed values for the radius, modulo the system in- 
clination, to explain the 1.46% amplitude ellipsoidal variations in 
J1741 and the 1.56% amplitude ellipsoidal variations in J2119. Our 
estimates include different values of the limb-darkening coefficient, 
but we assume a purely radiative WD atmosphere (ri = 1.0). For 
the 0.16 Mq elm WE) in J1741, its companion would be a neu- 
tron star if the inclination is less than about 64°, as marked by the 
vertical line, and if the inclination is less than 41°, the companion 
would exceed 3.2 Mq. For the 0.17 Mq ELM WD in J2119, its 
companion would be a neutron star if the inclination is less than 
about 49°. 

mary as a function of the orbital inclination of the sys- 
tem. Here we assume a linear limb darkening law and set 
ui = 0.36, but consider values between ui = 0.1 — 0.5, 
given the observations of [Parsons et al.l ([201 1[ ). Since we 
are most interested in setting a lower limit on the radius 
of the ELM WD, we also assume that the stellar ffux is 
purely radiative, and set ti = 1.0. This assumption may 
be well justified, as all models for 10,000 K low-mass 
WDs down to 0.17 Mq have no surface convection (J. 
Panci 2011, private communication). 

Figure [5] shows the radius of J 1741 as a function of the 
unknown inclination angle. Here we calculate the radius 
using three different values of the limb- darkening coeffi- 
cient. Our photometry rules out 2% or deeper primary 
eclipses, which constrains the system inclination to be 
less than 83°. The minimum radius of the ELM WD in 
J1741 is thus 0.074 Rq. This also implies that the mass 
of the unseen companion exceeds 1.11 Mq. 

The temperature and surface gravity estimates for 
J 1741 do n o t ove rlap with the theoretical models of 
[Panel et all ([2007[). making our radius measurement 
more important. I'Panei et al.[ ([2007| ) predict a radius of 
0.09 Rq for a 1 Myr old 0.16 Mq WD with Teff = 9,200 
K and \ogg = 5.7 cm s^^. Our radius measurement for 
J1741 provides the first observational constraints on such 
a low surface gravity WD. 

Using the values that best reproduce the cos2 ( ji term 
for ellipsoidal variations from [Morris fc Naftilad ([1993| ). 
the expected amplitude for the cos (j) term is less than 
0.02%. This is not compatible with our observations: we 
measure 0.23 ±0.06% variations in the cos^ term, which 
is a 3.8cr detection. 

This variation indicates that the primary is brighter 
by 0.23% when its heated surface is facing us at (/) = 0. 
Assuming that the unseen secondary is a 1.2 MqC/0 
WD with a typical radius of 0.0087 Rq {i = 75°), the 
secondary needs to be hotter than 78,00(3 K to reproduce 
a 0.23% reflection effect (,KopaL1959i) . Such a companion 
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Fig. 3. — High-speed photometry of J2119 over 10.4 hr (top 
panel), folded at the orbital period, 124.95 min. The bottom panel 
shows the same light curve binned by 100 phase bins. The solid 
line showrs a 1.51% amplitude signal at half the orbital period plus 
a 0.20% Doppler beaming effect at the orbital period. 

would easily be detected in the SDSS photometry and the 
MMT spectroscopy data. Hence, the companion may be 
a pulsar. 

Millisecond pulsars are known to "zap" their compan- 
ions, e.g., the b l ack w idow pulsar system PSR 1957+20 
(|Fruchter et al.l I1988D and the double pulsar system 
J073 7-3039 (|Burgav et all [200l iLvutikov fc ThompsonI 
|2005( ). This irradiation, which can lead to ablation of the 
companion star, is believed to be caused by relativistic 
particles emitted by the pulsar. The observed reflection 
effect suggests a pulsar companion, which would require 
i < 65°. If the inclination of the system is roughly 60°, 
the observed ellipsoidal variations would constrain the 
radius to > 0.078 Rq. A pulsar companion is unlikely to 
have a mass exceeding 3.2 Mq, putting an upper limit 
on the ELM WD's radius of 0.102 Rq (if ui = 0.1). 
Although not suggested in the models, we may also con- 
sider the scenario where some of the ELM WD surface is 
convective (ri = 0.48), in which case gravity darkening 
will increase the maximum allowable radius to 0.112 Rq 

if 7712 = 3.2 Mq. 

3. SDSS J211921. 96-001825. 8 

The ELM WD J2119 has T^s = 10,360±230 K, l og q = 
5.36 ± 0.07, and M - 0.16 Mq ()Brown et al.|[2Q10l ). It is 
a binary with a 2.0825 ± 0.0010 hr orbital period and a 
mass function of 0.501 ± 0.016 Mq, yielding a minimum 
mass for the companion of 0.75 Mq. 

3.1. The Light Curve 

Observations of J2119 were carried out at the McDon- 
ald Observatory over three nights in 2011 July, for a total 
of 1,252 points over more than 10 hr. All observations 
of this g = 20.0 mag WD were with 30 s exposures, and 
were reduced in an identical manner as those of J 1741. 

Figure [3] shows the Argos light curve of J2119 (top 
panel) folded over the best-fit orbital period from the 



radial velocity variations, 124.95 minutes. The bottom 
panel shows the same light curve binned into 100 orbital 
phase bins. The red solid line shows our best-fit model, 
as discussed below. 

3.2. Analysis and Discussion 

We model the light curve of J2119 using a five param- 
eter model as in J1741. Given the estimate for J2119 
of Tcff = 10,360 ± 230 K and a radial velocity ampli- 
tude K = 383 ± 4 km s-\ we expect a 0.41 ± 0.02% 
Doppler beaming signal. Only the cos 20 term for ellip- 
soidal variations and the sin0 term for Doppler beam- 
ing appear with any statistical significance in the Monte 
Carlo analysis of the Argos light curve (our full results 
can be found in Table[T]). This analysis yields an estimate 
for the observed amplitude of the ellipsoidal variations 
at 1.56 ± 0.13% and that of Doppler beaming signal at 
0.20 ± 0.12%. Argos photometry rules out 3% or deeper 
primary eclipses, which constrains the system inclination 
to be less than 85°. 

Again, we have computed a Fourier transform of our 
time-series photometry as an additional check. The 
highest peak in our dataset occurs at 3,745.9 ± 2.7 s 
(1.47 ±0.14% amplitude), at half the orbital period. Af- 
ter pre-whitening by this peak, the second-highest alias 
of the remaining peak occurs at 7,592±25 s (0.65±0.14% 
amplitude) . 

The disparity between amplitudes for the Doppler 
beaming signal in our Monte Carlo and Fourier analy- 
sis is notable, but considering that our photometric data 
cover less than five orbits for this relatively faint system, 
we will not over-interpret our results for this signal. We 
require more data to use the signal as a constraint on the 
radial velocity amplitude and spectroscopic temperature. 

On the other hand, the ellipsoidal variations are de- 
tected with high significance in both analyses. Again, 
we use the ellipsoidal variations to constrain the radius 
modulo the inclination of the system. We use the same 
limb- and gravity-darkening coefficients as in J1741, and 
our constraints on the radius of the ELM WD in the 
J2119 system are presented in Figure [21 Assuming that 
the ELM WD is purely radiative, we find a minimum 
radius of the ELM WD primary to be 0.10 Rq. This ra- 
dius is slightly larger than the radius estimate for J 1741. 
J2119 is estimated to be 600 K warmer than J1741. Even 
though this is only a relatively small temperature dif- 
ference, O.I6M0 WDs a re predicted to co ol down only 
by 1900 K after 10 Gyr (jPanei et al.ll2007[ ) and decrease 
in size by half. Hence, the relatively small temperature 
difference between J2119 and J1741 may be enough to 
explain the difference in their inferred radii. 

4. CONCLUSIONS 

We present the discovery of two new tidally distorted 
ELM WDs. We use the amplitude of the observed ellip- 
soidal variations to constrain the radius of this rare, fiuffy 
class of WDs. J1741 and J2119 are larger than 0.074 and 
0.10 Rq, respectively. These two obj ects plus Keple r Hot 
White Dwarf 3 (i? - 0.15i?Q iCarter et all [20ll and 
GALEX J171708. 5-1-675712 (R ~ 0.1J?.T^ IVennes et aP 
2011) are the largest radii WDs ever observed. Unlike the 
typical Earth-sized 0.6 Mq C/O-core WDs, these --0.2 
Mq He-core WDs are similar in size to a giant planet like 
Jupiter. 
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Fig. 4. — Observed mass and radius determinations for low-mass (He-core) WDs. The radius measurements for points marked by squares 
came from eclipsing systems and the dots from the amplitud e of ellipsoidal var iations, as in this work. The points for J1741 and J2119 
represent the case where i = 60° . The mass-radius tracks from lPanei et all |2007) cover a range of temperatures, ranging from nearly 22,000 
K, as seen in CS S 41177A, to nearly 8 ,500 K, as seen in NLTT 11748. Our mass/ra dius values come from several re ferences: (1) NLTT 
11748 (Stcinfadt et al.ll2010l: IKilic et al.t.20ia 'l: (2) KOI81B (van Kcrkwiik ct al. 2010VRowe et al.ll2010l); (3) KOI74B Hvan Kerk wiik et all 
[20TOh : (4, not pictured; 0.26±0.04 Mq, 0.15 ±0.0 1 Rp,) KH WD3 (Carter ct al. 2011); (5) J0106 (Kilic ct al. 2011b); (6 and 7) CSS 41177 
A and B (P arsons et a l. 20f^; (8) GALEX J1717 (IVennes et a l. 20111; (9) J0651 (Brown c t al.„2011b. ); (10) KHWD4 (Breton ct al...2011i) ; 
and J1741 and J2119 (this work). We have included the radii of Earth and Jupiter for reference. 



We have detected a Doppler beaming signal in the 
photometry of both systems, as a resuh of the high- 
amphtude radial-velocity variations. The effect mea- 
sured in J1741, with 0.59 ± 0.06% amplitude, represents 
the strongest Doppler beaming signal detected so far in 
a binary system. 

J1741 remains an intriguing object. It is so far the 
only tidally distorted WD that also shows a small but 
significant refiection effect. We observe a 0.23 ± 0.06% 
flux modulation at the orbital period of this ELM WD 
which peaks at (/> = 0. No comparable signal is observed 
in J2119, which shows even higher-amplitude ellipsoidal 
variations than J1741. Follow-up photometry at a larger 
telescope is required to confirm this reflection effect. If 
real, this may be due to the irradiation of the WD by 
a pulsar companion that may be detected in the radio 
or X-ray observations. We are pursuing follow-up X-ray 
observations of J1741 with the XMM-Newton telescope. 
If the companion is a neutron star, the inclination of 
J1741 is less than 65° and the WD has a radius > 0.078 

Rq. 

We measure ellipsoidal variations in J1741 and J2119 
of 1.46 ±0.06% and 1.56 ±0.13% amplitude, respectively. 
A Fourier analysis of these high-amplitude variations do 
not indicate a detectable difference between the period of 



the radial velocity variations and twice the period of the 
ellipsoidal variations (see Table [T|), although there is too 
much aliasing in these observations to say with certainty 
that the system is synchronized, such that the ELM WD 
rotation period matches the orbital period found from 
the radial velocity variations. 

Using these ellipsoidal variations, we have put limits 
on the radii of two of the lowest-mass WDs ever directly 
observed. There are less than a dozen other empirical 
mass-radius determinations for low-mass (He-core) WDs, 
and we plot our results with the theo retical mass-radius 
relations from evolutionary models bv lPanei" eHIIlpOOTl ) 
in Figure m 

Evolution models predict that He-core WDs with 
masses less than 0.17 Mc^ should sustain sufficient hy- 
drogen shell burning (iSerenelli et al.l l2002t iPanei etaP 
120071 iSteinfadt et"al1 12010( 1. In fact, a majority of the 
flux from these < 0.17 Mq WDs comes from this resid- 
ual burning. This should in turn cause a considerable 
increase in their radius, which has so far been borne out 
in our empirical mass-radius determinations. 

The orbital periods in both J1741 and J2119 are de- 
caying due to the emission of gravitational radiation. It 
is possible to use the time-of-minimum of the ellipsoidal 
variations in both objects to detect this orbital period 
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TABLE 1 

Physical parameters 



J1741 


RA (J2000) 


17:41:40.49 


Dec (J2000) 


+65:26:38.7 


Teff (K) 


9,790 ± 240 


logg (cm s~^) 


5.19 ±0.06 


mi (Mq) 


0.16 ±0.01 


90 (mag) 


18.271 ± 0.022 


Port (s) 


5,279.90 ±0.86 


Pelv (s) 


2,639.920 ±0.016 


COs((/))amp (%) 


0.23 ±0.06 


COs(2()!.)amp (%) 


1.46 ±0.06 


Sin((^)amp (%) 


0.59 ±0.06 


Sin(20)amp (%) 


0.01 ±0.06 




0.074 — 0.112 


'Po,ELV (BJDtdb) 


2455686.76251(21) 



J2119 


RA (J2000) 


21:19:21.96 


Dec (J2000) 


-00:18:25.8 


Teff (K) 


10,360 ±230 


logg (cm s~^) 


5.36 ±0.07 


mi (Mq) 


0.17 ±0.01 


90 (mag) 


20.000 ±0.021 


Porb (s) 


7,496.9 ±3.5 


Pelv (s) 


3,745.9 ±2.7 


COs((/))amp (%) 


0.02 ±0.13 


COs(2(/.)amp (%) 


1.56 ±0.13 


Sin((/))amp (%) 


0.20 ±0.13 


Sin(2(/.)amp (%) 


0.06 ±0.13 


ri (Ro) 


0.10 — 0.13 


T'o,ELV (BJDtdb) 


2455769.84065(66) 



decay. This rate of the orbital period change is sensitive 
to the mass of the unseen secondary, but for J1741, the 
orbital period is changing by at least —3.6 x 10~^^ s s~^. 
This would cause the observed minimum to occur more 
than 8 s sooner within 15 years. With less than 10 hr 
of data spread over four months, the formal error on the 
observed phase is roughly 20 s, which will make this mea- 
surement difficult in less than a decade. If the companion 
is indeed a pulsar, the orbital period would change by at 
least —4.5 x 10~^^ s s'-^, and the timing measurements 
of the pulsar could make this measurement possible in 
just a few years. For J2119, the orbital period is chang- 



ing by at least —3.4 x 10"^'' s s~ . We have included 
our To ephemeris measurements for this first epoch of 
observations in Table [T] 



The observations of J2119 in 2011 July were assisted 
by G. Miller, S. Wang, G. Earle, J. Pelletier, and A. Rost, 
undergraduate students at the University of Texas Fresh- 
men Research Initiative. J.J.H., M.H.M. and D.E.W. 
gratefully acknowledge the support of the NSF under 
grant AST-0909107 and the Norman Hackerman Ad- 
vanced Research Program under grant 003658-0252-2009. 
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